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Abstract 



Using transfer lots, we can overlap the processing of a batch on several consecutive 
machines, and thus reduce the makespan considerably. This in turn promotes work-in- 
process reduction. In this paper we investigate the transfer lots sizing problem for a given 
batch size under two operating procedures. Our objective is to minimize the makespan 
subject to a transferring budget. An important part of the solution involves partitioning the 
problem to subsets of machines without losing optimality. For each part (subset), the first 
and the last machines operate continuously while intermediate machines may idle intermit- 
tently. The first operating procedure we consider calls for the lots to be identical across 
all machines in each subset. The second operating procedure allows sub-lots for some of 
the machines or for some of the lots. Tliough more elaborate, the second operating 
procedure yields demonstrably superior results. The techniques provide satisficing feasible 
solutions, which can also serve as efficient bounds for an exact branch and bound integer 
linear programming model. 
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Introduction 



In recent years the Japanese have achieved monumental industrial success by implementing the just-in- 
time (JIT) production system on a nationwide scale [9; 4, pp. 736-769]. A basic tenet of JIT is that large 
batches are contra-productive in more than one way. For instance, they cause excessive work-in-process 
(WIP), excessive lead-time, and reduced flexibility. Large batches also compromise quality, because by 
the time a defect is detected it is too late to do anything about it. Therefore, JIT calls for small batches, 
ideally of one unit each. 

Other important elements of JIT-beyond the scope of this paper--are total quality management, 
workers’ participation (Quality Circles), and striving for constant improvement. Our main concern here 
is with aspects of materials flow. 

Also known as The Toyota Method, JIT is designed primarily for the repetitive manufacturing 
environment. It is a pull system, that is, usage downstream authorizes fabrication upstream. Assembly 
lines, often found in the repetitive manufacturing environment, are conducive to moving parts one-by- 
one, as urged by JIT. Parts required for assembly or fabrication are fed to the right stations in small 
containers. The units in each container usually make up a production batch. To avoid disruptions, 
buffers comprising a small number of containers are allowed in front of all stations. To avoid excessive 
WIP, strict limits on the number of containers in each buffer are observed. Part of JIT is a continuous 
effort to reduce these buffers, and still maintain smooth output. 

In the mass production environment there are few potential setups for each machine (we use the 
term machine as a generic for any station where the products have to be processed). To make small 
batches possible, these setups have to be vigorously streamlined. Reducing setups that used to take 
several hours to less than 10 minutes is a must under JIT. For an illuminating text on this issue and its 
impact on the evolution of JIT, see Shingo [16]. 

In contrast, for medium volume production, and even more so in custom job shops, a large variety 
of products are produced. Therefore, the number of potential setups increases, and it becomes 
progressively uneconomical to reduce all of them. Under such circumstances, specifying large batches 
may be necessary. 

Can we capture the major advantages of JIT--such as reducing the lead-time and the WIP 
inventory-without setting the machinery up more than once per batch, while still specifying sizable 
batches? Goldratt, the developer of OPT (Optimized Production Technology) [7; 12, pp. 692-715; 10], 
answered this question in the affirmative. Although OPT does not live up to its name, it is a 
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sophisticated production control system that successfully applies many JIT ideas to batch production. 

It is possible to adopt the OPT philosophy, also known cis synchronized manufacturing [4, pp. 790- 
839], without using any computerized system. Nevertheless, many perceive OPT as a competitor of MRP 
(I and II) [13; 12, pp. 655-658]. Our stance in this paper, following Vollman [20], is that OPT is a 
potential enhancement to MRP. There are four key features in OPT that most MRP packages do not 
support [10; 20]: (i) concentrating on bottleneck resources; (ii) scheduling activities on bottlenecks (and 
downstream from bottlenecks) forward instead of backwards, thus utilizing them fully; (iii) specifying WIP 
buffers (only) in strategic locations (for example, in front of bottlenecks); and, (iv) allowing transfer lots 
to be smaller than the batches they belong to, thus overlapping the processing on sequential machines. 

It is the transfer lots scheme that yields the major lead-time and WIP reductions that OPT 
achieves. According to a broad interpretation of the OPT principles, these lots need not necessarily be 
of equal size. Judging by the output of OPT, however, it seems that they do use transfer lots of constant 
size [10]. (In this paper we allow the lots to vary.) 

Although our main concern here is with (iv) and not with (i) through (iii), we note that the 
literature on scheduling is oriented to forward-scheduling, so it applies to scheduling bottlenecks [e.g., 
3; 5]. Also note that linear programming can be used not only to identify bottlenecks (i.e., binding 
constraints), but also to optimize the product-mix. Ronen [14] gives an analytic model for (iii), based on 
the newsboy model. Ronen and Starr [15] discuss the relationship between OPT and well-known 
optimization methods. 

Some work has also been published in the realm of (iv). Recent examples are Graves and 
Kostreva [8], and Truscott [18; 19]. The interested reader may find references to earlier efforts there. 
The assumptions in [8] are: (i) constant demand, and (ii) equal production rates for all machines. The 
model in [8] is developed for two machines. It optimizes the number of lots under the constraint that 
they should be strictly equal and integral. If more than two machines are involved, the authors apply 
their model on a pair-by-pair basis. [19] is based on [18], and does not make the constant production- 
rate assumption. Several stages are investigated under a restriction that once a batch starts on a machine, 
it is run continuously to completion. Transfers are limited to multiples of equal-sized sub-batches. 
Limitations on the transportation capacity are also taken into account. [18] and [19] are oriented towards 
implementation and dwell less on theoretical issues. 

Trietsch [17] obtains optimal lots for one batch on two machines. The assumptions are that (i) 
the units can be transferred one-by-one or in any combination, up to and including the whole batch; (ii) 
the batch size is given; and (iii) the number of transfers is either constrained by a budget or by a 
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limitation on the transportation resources (e.g., there are j vehicles available for the transfers, so lot j + 1 
cannot be moved until the first vehicle returns). The solution is then extended to several batches that 
have to be processed on the same two machines in the same order. Finally, a fast heuristic is introduced 
to extend the model to several machines on a pair-by-pair basis (similarly to [8; 18; 19]). In the latter 
case, transferring a lot incurs a cost that may be different for different machines, and there is a budget 
constraint on the total transferring expenditure. The number and composition of lots may change for 
each machine, to utilize the budget better. 

Let us look at a simple example: we have to process 250 widgets on 4 machines. Machines 1, 2, 
3 and 4 take 1, 2, 1 and 3 minutes per widget, respectively. Without splitting the batch to lots, the make- 
span is 1750 minutes. Suppose now that the budget allows two transfer lots from each machine. If we 
stipulate equal lots and require each machine to process all the widgets continuously (as in [8], but note 
that the production rate is not equal), the makespan will be 1375 minutes, a reduction of 21.4% 
(Figure 1). 

Insert Figure 1 about here 

By allowing intermittent idling in Machine 3 (as OPT does), while still specifying equal transfer lots (as 

OPT probably does), we can reduce the makespan further to 1250, a total reduction of 28.6% (Figure 2). 

✓ 

Insert Figure 2 about here 

Instead of allowing intermittent idling, if we allow the lot sizes to vary for each machine, as in [17], we 
can reduce the makespan to 1230, a total reduction of 29.7% (Figure 3). (See Section 6 for 
computational details of this and the following illustrations.) 

Insert Figure 3 about here 

In this paper we generalize the results of [17] for one batch and several consecutive machines by 
allowing intermittent idling. As in [17], we also allow varying lot sizes. In the present example this can 
further reduce the makespan to 1155, a total reduction of 34% (Figure 4). 

Insert Figure 4 about here 

A crucial part of the solution involves schemes designed to partition the problem without losing 
optimality. For each part (subset), we constrain the first and the last machines to operate continuously 
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(see Section 2 for a formal definition of partitions). For instance, Machine 2 in Figure 4 is the last 
machine of the subset {1, 2}, and the first machine of the subset {2, 3, 4}. The first operating procedure 
we consider calls for the lots to maintain their composition across all machines in each subset, but allows 
different lots across subsets, as in Figures 3 and 4. The second operating procedure allows the use of 
sub-lots. That is, we distinguish between parent lots that remain intact aaoss all machines in each subset 
as before, and sub-lots that make up the parent lots. Figure 5(b) illustrates such a case, where the first 
parent lot is one unit, and the second parent lot includes four units. On Machine 1, the second parent 
lot includes two sub-lots; Machine 2 recombines them to a single lot. 

Insert Figure 5 about here 

Given a large enough budget, the first procedure can always achieve any feasible makespan. If 
necessary, we can do that by transferring the units one-by-one. The second procedure, however, may 
achieve the same makespan with a smaller budget than the first, and in this sense it is superior. 

This paper develops fast satificing solution algorithms for minimizing the makespan under a 
transferring budget. The algorithms are easy to program and fast; their worst case complexity is 
polynomial. They yield feasible integral solutions, and are intended to be implemented in new or existing 
MRP systems. The algorithms do not require the use of any external mathematical programming 
packages. 

It is also possible to find the minimal makespan by integer linear programming (ILP). When 
solving by ILP, the satisficing solutions can serve as efficient upper bounds. An ILP model is presented 
in the Appendix. 

Following presentation of an early version of this paper in ORSA/TIMS St. Louis (October 1987), 
this author became aware that Baker was independently developing a similar model [1]. Baker restricts 
the lots to retain their composition across all machines (in contrast to retaining their compositions in each 
subset only here). He solves for two machines and several lots, and for three machines and two lots. His 
model approximates the solution by relaxing the integrality constraints on the lot sizes. The solution is 
by a set of rules inspired by a linear programming (LP) formulation. (The same formulation can serve 
to solve for several machines and several lots). 

Baker and Pyke [2] extend the results of [1] to several machines and two transfer lots, under the 
same assumptions. The solution is achieved by minimizing the maximal path in a network. [2]’s result 
for the first example would be lots of 107.143 and 142.857. The makespan is 1178.571, or 23.571 more 
than in Figure 4 (24 when integrality constraints are introduced). 
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The rest of the paper includes 11 sections. Section 2 introduces the formal problem. Sections 3 
through 5 deal with the first operating procedure: Section 3 examines a basic model for the first operating 
procedure, under the assumption that all intermediary machines can handle the loads assigned to them 
without requiring a partition of the problem; Section 4 investigates when and how to partition the 
problem; and Section 5 finds the minimal number of lots necessary to achieve the minimal makespan. 
Section 6 gives some simple examples. Section 7 introduces an exponential model with sub-lots, and 
Section 8 develops polynomial heuristics for the same purpose. Sections 9 and 10 take care formally of 
the issues of setups and integrality respectively. Section 11 introduces modifications that may be required 
for applying the model in practice. Finally, Section 12 concludes the paper with a brief list of related 
research questions. 



2. The Formal Problem 

[P] A batch of m items has to be processed sequentially on n machines, Mj, M2, ..., Each item 
requires T* time units of processing on Mj; for all i. Prior to processing the first unit, M. requires a setup 
time of SUj; for all i. Transferring a lot of any size (up to and including m items) from Mj to j costs 
Cj, and takes TTj time units; i = 1, 2, ..., n-1. It is required to minimize the makespan subject to a 
budget constraint on the total transferring expenditure, B. ■ 

Definition; The symmetric problem is obtained from [P] by reversing the order of the machines. In this 
paper we use the term symmetry to refer to the relationship between [P] and the symmetric problem. ■ 

Setups in [P] become tear-downs in the symmetric problem; otherwise, the symmetry is perfect. 
The symmetric problem has the same minimal makespan as [P], and can be solved by the same lots--in 
reversed order. Symmetry is instrumental in proving most of the results below, starting with the next 
theorem: 

Theorem 1 : Any feasible makespan can be realized in such a manner that both Mj and will process 

the whole batch continuously, although intennediary machines may have to idle intentiittently. 

Proof; Trivial for Mj, and by symmetry for Mj^. ■■■ 

Borrowing PERT/CPM terminology, the theorem simply suggests adopting "early start" on M j and 
"late start" on M^^. 

An item for our purpose may actually be a set of several units, say a dozen, if the policy is to 
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produce and transfer in dozens. For convenience, assume that at time 0 all the machines are free, but 
not yet set up for the batch. This assumption is not restrictive: if M- is busy at time 0, say until time t, 
simply add t to SUj. Similarly to [19], we also assume that transporting the lots is done independently 
of operating the machines. That is, the machines can continue working while the lots are being 
transferred. This assumption is appropriate in environments where dedicated resources are assigned to 
transferring items between stations. It is also appropriate if the transfer time, TT, is negligible. In 
addition, it may be possible for the operators of machines that idle between lots to handle the transfers. 
Other assumptions about this issue exist in the literature; e.g., see [6]. 

We use the budget constraint as an approximate way to allocate transportation resources to the 
various machines. In an environment where many such transfers are called for, and transportation is 
handled by a central department, this is equivalent to treating the transportation department as a profit 
center that sells transportation services to the jobs. 

[17] includes an analytic model where a prespecified number of vehicles are available. The 
solution specifies lots that are large enough to make it possible for the vehicles to return in time for the 
next transfer. This solution can be easily implemented for adjacent pairs of machines in the present 
problem. Nevertheless, it is difficult to generalize this solution when the same vehicle can serve more 
than one pair. 

The major effect of TTj on the makespan is increasing it by a constant, namely iT l j. This is 
true since we do not specify that the same vehicle has to handle all the transfers. Therefore, there is no 
need to wait for a vehicle to return from its former transfer before dispatching the current lot. 

In addition, TTj may influence the issue of whether machines downstream can be set up in time 
to process the first item that reaches them. (If they don’t, a binding constraint is introduced.) Until 
Section 9 we assume SUj = 0; for all i; therefore, for determining the optimal lots, we can also assume 
TTj = 0; for all i. 

What is the potential for makespan reduction here? The minimal makespan can always be 
realized by transferring the items one-by-one. Therefore, the minimal makespan is the processing time 
on the slowest machine plus the processing time of one item on all the others. Now subtract this from 
mTTj to obtain the maximal makespan reduction (MMR): 

MMR = (m - l)(2Tj - Max{Tj}). (1) 

In the example illustrated in Figures 1 through 4, the MMR is (250 - 1)(1 + 2+ 1 + 3- 3) = 996. 
Therefore, the reduction in Figure 1, 375, is 37.7% of the MMR. Similarly, in Figures 2 through 4 the 
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reductions are 50.2%, 52.2% and 59.7% of the MMR respectively. The reduction in the example 
illustrated in Figure 5(b) is 100% of its MMR. 

Similarly to [18], our stance in this paper is that as long as the makespan is minimized it is 
preferable to use as few transfers as possible. In this spirit we will show that 0(log m) transfer lots will 
often suffice to achieve the maximal makespan reduction even if the budget is not binding. 

We conclude this section with a definition of partitions. 

Definition : A set of machines that are required to work continuously is called a partition set, or simply 
a partition. We stipulate that a partition set must include Mj and Mj^ (as per Theorem 1). A partition 
can be read as a list of machines, or as a list of pairs. For instance, if Mj and Mj are adjacent to each 
other in a partition we say the partition includes the pair (i, j). For convenience, we list a partition by 
only listing the indices of the machines in it. When a partition includes r machines in addition to Mj and 
Mj|, we may list it as {p(o) = l, p(l), p(2), ..., p(r), p(r+l)=n}, or {p(s)}s=o,r+l- " 



3. A Preliminary Model 

In this section we start treating a relaxed problem, where fractional items can be transferred. The relaxed 
problem is solved by the relaxed solution, as opposed to the integral problem/solution. To avoid an 
excessive gap between the relaxed solution and the integral one, we stipulate that in a relaxed solution 
all lots should be >. 1. This restriction is also instrumental for demonstrating that the number of lots 
required to achieve the maximal possible makespan reduction is often 0(log m). We refer to instances 
where lots are allowed to be < 1 as super-relaxed. 

Using Theorem 1, we specify that Mj and M^^ should operate continuously, while the intermediate 
machines may idle between lots. We use our first operating procedure, i.e., the lots retain their composi- 
tion across all machines unless a partition is involved. We denote the size of lot j by Lj; j = l,...,k; we 
may also use L. informally as the name of lot j. We denote the cumulative sum of the first j lots by S-; 

J J 

e.g., Sj = Lj, and Sj^ = m. Under partition, Lj may be different for each subset, so formally we should 
use a double index to identify Lj and Sj. Nevertheless, it is possible to use the simpler notation without 
causing confusion. Our formal problem is now: 

[Pk] Solve [P] under the following assumptions: (i) TT. = SUj = 0, for all i; (ii) the lots retain their 
composition across all machines in each subset; (iii) fractional items may be transferred, but Lj >. 1; for 
all j. ■ 
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Definition: A solution is called well-behaved if (i) Mj and operate continuously as per Theorem 1; 

and (ii) Mj (i = 2, n-1) can process each lot as soon as it becomes available from (i.e., lots are 

neither delayed at nor queued at M.). If a lot can reach a machine before the machine is ready for 
it, we say there is a conflict. ■ 

Figures 3 and 4 illustrate well-behaved solutions. Figures 1 and 2 illustrate conflicts; e.g., in both 
cases the second lot from Mj can reach M 2 at time 250, but M 2 is not ready for it until time 375. 

It turns out that if fractional items may be transferred, the optimal solution tends to be well- 
behaved. The only exception may occur due to the restriction Lj ^ 1, which constitutes an integrality 
constraint when it is binding. Thus, the optimal super-relaxed solution is well-behaved. 

This is true because if conflicts exist, lots which have to wait could have been increased without 
increasing the makespan. Since the sum of all the lots is constant (m), some preceding lot could have 
been decreased, thus feeding Mj^ sooner, and reducing the makespan. This argument fails if by reducing 
the preceding lot we violate the constraint Lj ^ 1 for it. 

We also assume in this section that partition will be neither specified nor required; i.e., there exists 
a well-behaved solution that can be found without specifying any partition. Without partition, under our 
operating procedure the maximal number of transfers allowed by the budget is k = INT(B/2Cj). To 
obtain a well-behaved solution we must have 



in a well-behaved solution. Figure 6 illustrates this point: Lj = 100 items starts processing on M 2 at time 
100 and finishes on M^ at time 700, a total of 100(1 + 2 + 3) = 600 time units; I -2 = 150 units starts 
processing on Mj at time 100 also, and finishes on M^ at time 700--just when M^ can start processing 
it--also a total of 150(1 + 1 + 2) = 600 time units. 




( 2 ) 



That is, the time it takes to process Lj + j on the first n-1 machines should equal the time it takes to 
process Lj on the last n-1 machines. Thus Lj reach M^^ exactly when it is ready for it, as it should 



Insert Figure 6 about here 



Define 



Qy - (Ti+, + +...+ Tj)/(Tj + Ti^, +...+ Tj.,), 

and let Q = Qj j,; then the set (3) replaces the set (2) 

Lj = Lj jQ = ; for all j ^ 2. 



( 3 ) 
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In addition we have 



Lj + 1-2 +...+ Lj^ = m. (4) 

In order to satisfy (4), we use the familiar geometric progression to obtain 
m(l - Q)/(l - Q*') ; Q^l 

U = { (5) 

m/k ; Q = 1. 

Along with (3), this suffices to specify all the lots. To continue, let us investigate the conditions under 
which (3) and (5) yield a well-behaved solution. The question is whether there will be conflicts at the 
intermediary machines. Checking the following n - 2 feasibility conditions gives the answer. 

^l,n ^ Qfyp ; 2 j< p j< n-1. (6) 

The conditions are equivalent to comparing the lot sizes obtained in the solution to those obtained 
by stipulating that Mp will work continuously, given Lj. Under this stipulation (3) determines L3, 
etc., but with Q = Qi „ instead of Q = _. As long as the lots are not smaller than those required 

by Mp to operate continuously, they cannot reach Mp before it is ready for them. Note that the 
conditions are invariant under the number of lots, the size of the first lot, and m. In the example 
illustrated in Figure 6, Qi 4 = Qi 3 = 1-5- That is, one of the feasibility conditions is satisfied as an 
equality; therefore M^ operates continuously even though no partition is enforced. 

If conflicts do occur in intermediary machines, we can specify that the machines involved shall 
operate continuously, as do Mj and Mj^. That is, we partition the problem (see Section 4). 

Let [(i)] denote the numerical value obtained by applying (i). If [(5)] < 1 or [(5)]Q*^‘^ < 1, we 
use the following procedure: 

The Adjustment Procedure; 

(a) Q ^ 1 AND [(5)] < 1: 

1. Set Lj = Sj = 1, and let i = 1 (recall Sj = Sj^jjLj); 

2. if S- < m, let = LjQ (as per (3)), Sj^j = min(S. + L-^j, m}, i = i + 1; proceed to the 
next step; 

3. if Sj = m, set K = i, Lj^ = m - and proceed to the next step; otherwise, return to the 

former step. 

4. iterating backwards from] = K, if Lj < 1 increase Lj to 1 by decreasing Lj 1; upon encountering 
j such that Lj ^ 1, STOP. 
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(b) Q < 1 AND [(5)]Q*‘*^ < 1: 

Apply Part (a) of the procedure to the symmetric problem. ■ 



If Q ^ 1, the lots under (3) (except possibly the last one when Lj f [(5)]) are non-decreasing, so 
if Lj = [(5)] ^ 1 all the lots are ^ 1, and there is no need for adjustment. A symmetric observation 
holds if Q < 1 and Lj = [(5)]Q*^'^ >. 1. Othenvise, it is easy to verify that the procedure ensures 
Lj .> 1; for all j. Let * superscripts denote optimal values, then 

Theorem 2; If the feasibility conditions are satiated, [(5)] ^ 1 AND [(5)]Q*^*^ >. 1, then L^ = [(5)]; else, 
if the feasibility conditions are satisfied, then The Adjustment Procedure yields an optimal solution. (See 
Appendix for proof.) ■■■ 



r ^ 

Let MS denote the relaxed makespan, and MS denote the optimal integral solution, then 
MS - MS = h ^ 0 is the difference in makespan due to the integrality constraints. We conclude this 
section by developing a simple upper bound for h*, still under the assumption that no partition is 
involved, and that SUj = 0; for all i. Under partition, the bound will apply to each part separately. The 
sum of the bounds of the parts will then be the bound we seek. To continue, let fj be the fractional part 
of Sj, and define 



1 - fi ; fj > 0 

0 ; fj = 0. 



Theorem 3: h 



<_ min{Max{ej}2. 

i * ‘ 



i,„.,Tj, Max{fj}2j,2,„Tj} 



< "i-.n'Tj 



- Ma|c{Tj, T„}. 



Proof; We introduce two heuristics: 

Heuristic 1; Round all the S- values up to the nearest integer, i.e., Sj + ej. That is, 
the new value for Lj will be (Sj + ej) - (Sj_j + Cj_j). 

Heuristic 2: Truncate all the Sj values down to the nearest integer, i.e., Sj - fj. 

See Appendix for proof that Heuristic 1 increases the makespan of the relaxed solution by 
Max{ej}2j^j„.jTj, Heuristic 2 by Max{fj}2j^2,nTj, and both yield feasible solutions. ■■■ 

Users who are satisfied with approximate solutions can use one of these heuristics instead of 
applying the results of Section 10. 
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4. 



Partitioning the Problem 



The simplest way to enforce the feasibility conditions (6), when some of them are violated, is to look for 
the largest index p which maximizes Qj p, and specify the partition {1, p, n}. This creates n - p - 1 new 
feasibility conditions downstream from Mp, i.e., 

Op^>Qp^;p+i<i<ri-h 

and more machines may have to be added to the partition set between p and n. In this section we 
investigate this and other partitions of the problem. 

Definition : A pair (i, j) is called feasible if Q: t >: Qj for all i ^ p ^ j. That is, M; can feed M; as 
per (3) and (5), with both of them operating continuously. ■ 

Note that adjacent pairs, e.g., (i, i+1), are always feasible. 

Definition : A partition is called feasible if all its pairs are feasible. ■ 

Specifying any feasible partition and using (3) and (5) for each part yields a well-behaved solution. 

Definition : The partition implied by the feasibility conditions is called the minimal partition. That is, 
if {p(s)}s=o,r+i denotes the minimal partition, then 

^P(s)>P(s + l) “ ^^j>p(s)^^p(s)J^’ all 0 ^ s ^ r, and p(s+l) is the largest index j that maximizes 
Qp(j)j. We refer to the minimal partition as minpartit, and note that it is feasible by construction. ■ 

Definition : The (feeisible) partition which includes all the machines is the maximal partition, denoted 
by MAXPARTIT. ■ 

Definition : The medium partition, is the feasible partition with the maximal number of machines in it, 
such that given a large enough budget it will be possible to realize MMR, as per (1). We refer to it as 
MEDPARTIT. ■ 

It may happen that specifying MEDPARTIT will cause the makespan to increase relative to 
MINPARTIT. The definition merely says that given a large enough budget MMR can be achieved. The 
first example in Section 6 illustrates this point. 

MAXPARTIT contains all other partitions, feasible or infeasible, minpartit, as the name suggests, 
is contained in any feasible partition (see Theorem 6). And, any feasible partition properly contained in 
MEDPARTIT enables MMR. We show how to construct MEDPARTIT and substantiate these claims below. 

To continue our investigation of partitions, we present a few theorems. The proofs are deferred 
to the Appendix. Theorem 4 is the motivation for the definition of MEDPARTIT. It indicates that 
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partitions may limit the maximal possible gain. This is not surprising, since they impose additional 
constraints on the problem. Before stating the theorem we define MMRy as the maximal makespan 
reduction for pair (i, j). As in (1), we obtain 

MMRy = (m - - Max{Ts},=ij). (7) 

Theorem 4 : For any partition {p(o) = l, p(l), p(2), ..., p(r), p(r+l)=n} (feasible or not), 

MMRj n > 2i,o,rMMRp(i) p(j^.i). ... 

It is easy to construct examples where the inequality is strict (see Figure 7). We also want to show 
that the minimal partition does not imply any potential losses in MMR. Thus, for the minimal partition 
Theorem 4 is satisfied as an equality. We state this result below eis Theorem 5. But first, we present a 
lemma. 

Insert Figure 7 about here 

Lemma 1 : Let and M^be any set of three machines such that i < p < jandTp > Max{Tj, Tj}, 

then (i, j) is an infeasible pair. ... 

A direct corollary of Lemma 1 is that at least one of the slowest machines must belong to any 
feasible partition, including MINPaRTit and MEDPaRTIT. 

Theorem 5 : Let MINPARTIT = (p(o) = l, p(l), p(2), ..., p(r), p(r+l) = n}, then 

MMRin = 2i=o,rMMRp(i)p(j^.j). ... 

Theorem 6 : Let PARTIT be any feasible partition which is not identical to MINPARTIT, then MINPARTIT must 
be properly contained in PARTIT. ... 

We also need another lemma, which specifies sufficient (but not necessary) conditions for a pair 
to be feasible. 

Lemma 2 : Let (i, j) be a pair such that Tj^ ^ min{T., Tj} ; for all i < k < j, then (i, j) is a feasible pair. 

... 

We are now ready to identify MEDPARTIT. The procedure is based on the results above, and 
observation of the proof of Theorem 4 (see Appendix). We are looking for the largest set of indices, {s}, 
such that (i) T^ is monotone non-decreasing for any s which is smaller than the largest index of a slowest 
machine; and (ii) T^ is monotone non-increasing for any s which is larger than the smallest index of a 
slowest machine. 
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In more detail, by Theorem 5 we can start with the minimal partition, and add machines between 
its pairs where possible. It is enough to concentrate on one pair at a time, say (p(s), p(s+ 1)). There are 
up to three types of pairs in MINPARTIT; (i) pairs where Tp(s) < Tp(s+i); (ii) pairs where 
'Tp(s) = which can only happen if both of these values are maxima; and (iii) pairs where 

^p(s) ^ ^p(s+l)- If we have a pair of the first or second type, we look for the smallest index i such that 
p(s) < i < p(s+l) and — ^p(s)> and include it in MEDPARTIT. By Lemma 2, (p(s), i) is a feasible pair. 

— ^p(s + l) (otherwise, by Lemma 1, the pair (p(s), p(s+l)) would not be feasible, and it must be 
feasible because it belongs to minpartit), so (i, p(s+ 1)) is also a pair of the first or second type (though 
not necessarily feasible). Thus we can repeat the procedure iteratively for the pair (i, p(s+ 1)), where i 
takes the place of p(s) above. When no such i exists, we move to the next pair of MINPARTIT. Pairs of 
the third type are treated symmetrically to pairs of the first type. 

Any feasible partition contained in MEDPARTIT also has the same structure, i.e., (i) T^ is monotone 
non-decreasing for any s which is smaller than the largest index of a slowest machine; and (ii) T^ is 
monotone non-increasing for any s which is larger than the smallest index of a slowest machine. 
Therefore, such partitions also enable MMR. 



Makespan Reduction as a Function of the Number of Transfers 



Let (i, j) be a feasible pair, specified as part of a partition, and let MRjj(k) denote the marginal 
makespan reduction associated with allowing k transfer lots relative to k-1 lots for the subset Mj,...,Mj. 
Let Kj • be the number of lots required to achieve MMR. • (as indicated by running The Adjustment 
Procedure for the subset, or see (11) and (12) below), then 
MRjj(l) = 0 (this transfer is essential). 



MRjj(k) = { 






; Q,j = 1, 1 < k < K,j, 



( 8 ) 



MRy(Kij) = MMRy - MRy(Kjj - 1). (9) 

(8) is derived from (A9) (see Appendix). Next, we have a convexity property: 



Theorem 7 : 77je {MR. .(k )}.^23 „ ''lonotone Jecreon'/ig. ■■■ 

Furthermore, the first extra transfer can yield at least half of the maximal makespan reduction, 
MMR; ;, i.e.. 
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Theorem 8 : MR; s(2) > MMR; ;/2. ■■■ 

In case of partition, Theorem 8 applies to each part separately, but not necessarily to the whole 
set. For instance, in the example illustrated in Figure 7(a), where MAXPARTIT is specified, the total make- 
span reduction is 18— less than half of MMR = 99. Nevertheless, it is easy to verify that under 
MEDPARTIT or any feasible partition contained in it (e.g., minpartit). Theorem 8 does apply to the whole 
(Figure 4). 



Identifying the Optimal Partition 

Given a feasible partition (e.g., minpartit, medpartit, or maxpartit). Theorem 7 suggests allocating 
the budget to the parts of the problem by a heuristic procedure, to which we refer as The Greedy 
Heuristic. First, we list and sort all the MRjj(k)/2j,jj.;Cg values. Then, we specify as many transfers 
as possible from the top of the list. That is, we use transfers that yield the highest marginal gain per 
dollar-as long as the budget allows. This also makes it possible to find the most economical budget, if 
the budget is not the result of real technical constraints. To that end, we simply add transfers from the 
top of the list as long as the extra expense is justified by the marginal contribution. If the heuristic 
allocates the budget completely, the solution is optimal. Otherwise, it may happen that by giving up some 
of the transfers indicated by the heuristic, larger makespan reductions may be obtained in other pairs. 
This is the heuristic presented in [17] for MAXPARTIT. The main problem with The Greedy Heuristic is 
that we need to know the partition in advance, or check all the possible candidate partitions (which would 
lead to an algorithm exponential in n). 

It is possible, however, to solve this problem to optimality by dynamic programming (DP). The 

Appendbc presents a DP algorithm that finds the optimal partition and allocation in 0(n m ) or 
'1 

0(n mlog m), depending respectively on whether or not Q; •, is likely to be 1. 



5. Achieving the Maximal Gain with the Minimal Number of Lots 

Let MSy be the minimal makespan possible for a feasible pair (i, j) when k transfers are allowed. MSy 
is measured from the start of the batch on Mj to its finish on Mj. Observe that we have to process the 
first lot on the first j - i machines and then process the whole batch on Mj. This leads to the following 
bound 



MS; 



y > Lj(Tj + Ti ^.1 +.... + Tj.i) + mTj (^ (T- + +.... + Tj.i) + mTj). 



( 10 ) 
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The second inequality arises because Lj >_ 1. By symmetry we also have 

MSjj ^ + T-^2 +--+Tj) + mTj ( ^ + T-^2 +— -+Tj) + mTj). 

Either or both of these inequalities can be strict, but under (3) and (5) they are both satisfied as 
equalities. For that to happen, both Mj and Mj have to operate continuously. Our objective in this 
section is to find the minimal number of transfers required to minimize MS- ;. We denote this number 
by Kjj (corresponding to K in The Adjustment Procedure). This number may be slightly impacted by 
the integrality constraints that we relaxed for the present, but it gives a good approximation for the 
integral version as well. Our procedure for finding K. . is simple. First, if Q: j > 1, then (3) leads to 
increasing lots, and Lj is the smallest lot. Since Lj ^ 1, by specifying Lj = 1 and using (3) for the other 
Lj values we obtain an increasing geometric series. (10) assures us that this policy will minimize the 
makespan, since it will hold as an equality. How many members should the series have to exceed m for 
the first time? The answer is: 



Kjj = SUPINT(log[(Qy - l)m + IJ/logQjj) ; Qy > 1, (11) 

where SUPINT(x) = smallest integer ^ x. Next, if Q: j < 1, then by symmetry 
K.. = SUPINT(log[(l/Qij - l)m + l]/log[l/Qij]) ; < 1. (12) 

Finally if Q: s = 1, then K. . = m. 

The procedure of calculating Lj by (3) based on L, = 1 is similar to our approach in The 

J 

Adjustment Procedure. Note that as a result, we obtain 
Lkj . = m - Sj^. J ^ 

The optimal solution is not necessarily unique here, since we have some degree of freedom in changing 
the lots following Lj without violating (2). Nevertheless, the lots specified above are optimal. We say 
this not only in the sense that the makespan is minimized but also in the sense that only necessary 
transfers are specified. Similarly, if the bound is on the last lot instead of the first, we use the same 
procedure, but backwards. Here we use (12) below instead of (3), starting with Lj^^. 

Lj(l/0) ; 3ij<Kij 

M-1 ~ ^ 

^ ni - ; j = 2. 

By observation of (11) or (12), if Q ^ 1, K- ; is 0(log m). This is obviously good news to 

‘J 

management. It also has a beneficial effect on the complexity of the procedures we employ. 
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6. Examples 



To illustrate the first operating procedure, along with its shortcomings, let us look at some simple 
examples. In Examples 1 and 2, with four machines each, we look for the minimal makespan with two 
transfer lots> Examples 3 and 4, with three machines each, illustrate the budget allocation and its effect 
on the makespan. 

Example 1 : Let m = 250, Tj - "^2 ~ ^3 ” "^4 “ assume the budget suffices for exactly two 

transfer lots. (We ignore the issue of optimal budget allocation here.) 

Solution (see Figure 6 ): Qj 4 = (T 4 + T 3 + T 2 )/(T 3 + T 2 + Tj) = (3 + 2 + l)/(2 + 1 + 1) = 1.5, 
Qj^ = (T 3 + T 2 )/(T 2 + Tj) = (2 + 1)/(1 + 1) = 1.5 ^ 1.5, and = 1 ^ 1.5, = = > no partition 
is required. That is, MINPARTIT = {1, 4}; MEDPARTIT = MAXPARTIT = {1, 2, 3, 4}. Since = Qj 4 , 
M 3 will operate continuously under MINPARTIT, even though it is not specified explicitly in the partition. 
Solving for Lj we obtain 100 units, and the total makespan is 1150. ■■ 

This example serves to illustrate an important point about MEDPARTIT; given a large enough budget, 
MEDPARTIT enables the maximal possible gain; but if we restrict the budget, as implied in this example, 
then partitioning at a machine which belongs to MEDPARTIT can increaise the makespan. In this example, 
there would be an increase of about 8 time units associated with adding M 2 to the partition. We leave 
the computational details to the interested readers. 

Insert Figure 8 about here 

Example 2 : As before, but interchange the order of M 2 and M 3 , i.e., T 2 = 2 and T 3 = 1. Here, if we 
try to stick to the former solution (Lj = 100), there will be a conflict at M 2 (see Figure 8 ), since M 2 still 
processes Lj when 1^ reaches it. Indeed, = 2 > 1.5 = Qj 4. The total makespan is 1200. By 
specifying MINPARTIT-(1, 2, 4} in this example--we can reduce the makespan to 1155 (see Figure 4), most 
of the way back to 1150. Note that this requires changing the size of Lj from 83 between Mj and M 2 
to 107 downstream. Figure 3 illustrates the makespan for the same example under MAXPARTIT, again with 
a variable size for Lj. The additional delay is 75 units, or about 6%. (For this example, the integrality 
constraints were satisfied by simple rounding, and therefore the makespan may be slightly off.) ■■ 

Example 3 (see Figure 7): Let m = 10, Tj = T 3 = 10, T 2 = 1, Cj = C 2 ® = 20. 

Solution Under First Operating Procedure : Each compound transfer costs Cj + C 2 = 2, so we can have 
up to 10 lots. Q = (1 + 10)/(10 + 1 ) = 1 . = = > Lj = 1 , and we obtain 10 equal lots of one item each. 
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Checking the feasibility conditions (6), 1 > 0.1, hence no conflicts will occur. The makespan is 111. ■■ 

Solution Subject to All Machines Operating Continuously : This is a forced partition, namely MAXPARTIT. 
Applying (12) to pair (1, 2), we obtain SUPINT(log(10 - 1)10 + 1]/Iogl0) = SUPINT(log91/loglO) = 
2. The same numeric result applies to (2, 3), due to symmetry. Hence, the budget is not binding. By 
(5) we obtain Lj for pair (1, 2) = 9/.99, but since this leads to 1-2 < 1, The Adjustment Procedure yields 
1-2 = 1, and Lj = 9. For the second pair, by symmetry, Lj = 1, and 1-2 = 9. The total makespan is 192 
(and not 191), since M 2 is constrained to work continuously, and hence it cannot start until time 91. M 2 
feeds Mj for the first time at time 92. M^, though it can start immediately, requires 100 time units. In 
this example Theorem 4 is satisfied as a strict inequality. Indeed note that M 2 does not belong to 
MEDPARTIT. ■■ 

Example 4 : Let m = 5, Tj = T 2 = 1, T^ = 5, Cj = 1, C 2 = 10, and B = 23. 

Solution Under First Operating Procedure (see Figure 5[a]): Q = (1 + 5)/(l + 1) = 3. Cj + C 2 = 11, 
so the budget suffices for two composite transfers, leaving $1 unused. By (5) Lj = 1.25, and the 
completion time is 27.5. Now, if we consider the integrality constraints, we are indifferent between 
Lj = 1 or Lj = 2. In the former case the second lot of 4 items will get to M^ at time 9, and it will still 
require 5*4 = 20 time units to finish the batch. In the latter case, Mj receives the first lot at time 4, and 
finishes at time 29 again. ■ ■ 

Solution Subject to the Maximal Partition : This is a forced partition again, but note that this time 
MEDPARTIT = MAXPAJtTlT. Qj^ = 1, SO we can use up to 5 lots there. Applying (11) to pair (2, 3), we 
obtain SUPINT(log(5 - 1)5 + l]/log5) = SUPINT(log21/log5) = 2. That is, there is room for up to one 
extra transfer there. Here our budget is binding. By (8) we get MR2^(2) = 4; MRj^(2) = 2.5; 
MRi^(3) - 0.8333; MRj^(4) = 0.4167; MRj^(5) = 0.25. Dividing these values by the respective 
transfer costs, the contribution per dollar of MR2^(2) is 0.4, while the other values remain unchanged. 
Therefore, the greedy heuristic will specify three extra transfers between 1 and 2. Then, it will try to 
introduce an extra transfer for (2, 3), but the remaining budget of $9 will be $1 short. In this example 
an exact algorithm (such as the DP model) can yield the optimal allocation, namely, one extra transfer 
for (2, 3) and two for (1, 2). The resulting makespan is 27.667. The best integral solution here yields 
28. (Note 27.667 > 27.5 but 28 < 29. Since 27.667 > 27.5, the dynamic programming algorithm would 
not specify a partition, even though with the integrality constraints it is better.) ■ ■ 

Optimal Solution (See Figure 5[b]); Allocate $3 to three lots from Mj, and $20 to two transfers from 
M 2 , thus utilizing the budget fully. Let Mj send the first item as the first lot, and then two lots of two 
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items each. M 2 will be able to send the first item to M 3 at time 2 , and start the next four items at time 
3 (M 2 will idle between 2 and 3). At time 7 M 2 and M 3 both finish their current load, and M 2 sends 
the four items to M 3 "just in time." Thus = 27. ■■ 



The last solution happens to comply with the integrality constraints and achieves the maximal 
makespan reduction as per ( 1 ). Therefore, it is the optimal solution. Note that this integral solution gives 
a shorter makespan than the 27.5 we got before with the relaxation. It compares even more favorably 
with the integral solution’s 29. Thus, our first operating procedure leaves something to be desired. 
Indeed it is intuitively clear that when the transfer costs and the potential marginal gain differ from stage 
to stage, forcing the model to specify the same number of transfers across all stages may be wasteful. 
This motivates our second operating procedure. 



7. Introducing Sub-lots into the Solution 

In this section we treat [Pk] without the constraints Lj ^ 1; leaving them in would complicate the 
presentation considerably. Thus, we look for the optimal super-relaxed solution. Section 10 develops 
integral solutions for this section’s operating procedure (as well as for the former operating procedure). 

Examine the optimal solution of Example 4. At first, the lot size remains constant for all machines 
(one unit). Then, two lots from Mj to M 2 are combined to one lot from M 2 to M 3 . Note that though 
M 2 has to idle after the first lot, it does so only after having dispatched all the items it processed. One 
way to describe the optimal solution of Example 4 is to say that the batch is divided to two parent lots. 
The first is processed on all machines "as is," while the second is further divided to two sub-lots on Mj. 
Each machine processes each parent lot continuously. The intermediary machine idles between parent 
lots. Note that the solution is well-behaved. 

This leads us to an operating procedure where we retain the same parent lots across all machines, 
but allow sub-lots within each parent lot. Machines can idle between parent lots, but not between sub- 
lots. Since the number of sub-lots in parent lot i and parent lot j need not be equal. Therefore the 
budget should not necessarily be allocated to the parent lots equally. Thus we are presented with a sub- 
problem of allocating the budget to the parent lots. We need to solve the sub-problem before allocating 
the sub-budget within the parent lot. Initially, we assume all intermediary machines idle between parent 
lots. By Theorem 1 we know that Mj and Mj^ do not have to idle. Later we will discuss partitions, and 
thus take care of conflicts in intermediary machines. 
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Assume for a while that we know how many parent lots there are, and to how many sub-lots each 
parent lot is divided on each machine. We refer to this information as the grouping information. We 
proceed to deduce the optimal lot (and sub-lot) sizes from the grouping information. 

The size of parent lot i, still denoted by L-, is constant across all the machines. The sub-lots may 
be different for different machines, even if their number is the same. Observe that under our assumptions 
each parent lot is processed as a batch under maxpartit. When a parent lot is not divided to sub-lots 
on some machine we may refer to it as a sub-lot of itself. Let: 



• k relative size of sub-lot j, emanating from M|^, and belonging to Lj. (^jYy = 1.) 

• ST- be the time the first sub-lot of L| starts processing on Mj^. 

Then, since Mj and operate continuously, the first sub-lot of L| starts processing on Mj at 

time 



and on Mjj at time 



(13) 



ST. 




Let Lj denote a tentative value for Lj, then to calculate all the L- and Y- ; 
1 and proceed recursively, as outlined below. 



(14) 

values we start with 



(1) At stage i, given l| and Y_ • for p < i and all j, k, calculate Y; • ^ for all j, k; if i = K, go to (3); 
else, set i = i + 1 and go to (2); 

(2) Given L* and Y„ j u for p < i and all j, k, calculate l| and return to (1) 

(3) For i = 1 to K let Lj = mL|/2L* 



The details of the calculations are deferred to the Appendix. By observing these details we can 
see that unless there are conflicts between parent lots on intermediary machines we obtain well-behaved 
solutions. 

We still operate under three assumptions: (i) we have the grouping information; (ii) all 

intermediary machines are allowed to idle intermittently; and, (iii) there are no conflicts at intermediary 
machines. Given these assumptions, we have an efficient procedure to calculate all the lot and sub-lot 
sizes. It is when lifting the assumptions that our method may become exponential. 

For instance, to lift the first two assumptions we can generate all the possible groupings and 
partitions which the budget allows. It is easy to see that this can be done in exponential time. We then 
solve for each such grouping/partition as described above. 
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As for conflicts at intermediary machines, we check for them by the following set of conditions; 



STi,k ^ STj.j + Lj.iT^ ; k = 2,..., n-1, all i. (15) 

Where we calculate the ST values recursively (using (13) for k = 2), i,e., 

^ ^ 

If the conditions hold, the solution is a valid candidate. Otherwise, that particular 
grouping/partition is rejected. The best candidate solution is the method’s final output. 



8. Polynomial Heuristics with Sub-lots 

In this section we introduce simple rules based on results from the former analysis, to obtain polynomial 
heuristics. 

For a while, let us assume that no partition is required. The technique we develop under this 
assumption can then be applied to the separate parts of the problem, if we partition it later. Further- 
more, it will indicate one possible "quick and dirty" method of partition. 

Under this assumption, we have to allocate the budget to parent lots, and-within each parent lot- 
-to the different machines. This gives us the grouping information. Given the grouping information. 
Section 7 lists all the necessary calculations. 

The basic rule we use for allocating the budget among the parent lots is to make the allocation 
as equal as possible. To the extent it is not possible to allocate the budget evenly, we make the allocation 
monotone non-decreasing, or monotone non-increasing. To choose between those two, we have the 
following rule of thumb: if the first machine is faster than the last, choose the non-decreasing order; 
otherwise, choose the non-increasing order. (This rule works well for Example 2. The rationale behind 
it is that it allocates more to the larger parent lots.) 

Allocating the budget evenly does not yet tell us how many parent lots we should use. Dividing 
the budget by 2Cj will give us an upper bound for that purpose. Another plausible (heuristic) upper 
bound may be obtained by using Ky, as per (11) or (12). Next we may check for lower numbers of 
parent lots. We continue considering less and less parent lots until one of the following happens: (i) we 
get down to a single parent lot; (ii) our results start to deteriorate instead of improving. Symmetrically, 
we can start with one parent lot, and then increase the number. 
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At this stage, it remains to allocate the budget within the parent lot. To do that we propose to 
use the solution of the maximal partition procedure as applied to the parent lot (instead of to the whole 
batch). Note that we do not really have to know m, or Lj, to use this solution. Hence we have enough 
data to apply it. Also recall that the partition is known (i.e., MAXPARTIT), and The Greedy Heuristic as 
well as the DP solution can supply the necessary allocation. 

Theoretically, by applying the maximal partition procedure we minimize the makespan of the 
parent lot, s.t. no idling between sub-lots. In fact, what we wish to do is to maximize the size of the 
particular parent lot, without exceeding the time frame allotted to it. That is, we wish to maximize the 
throughput. Minimizing the makespan for a known number of items, however, is dual to maximizing the 
throughput during a set period. Therefore, this solution is appropriate here. 

Having made the allocation to parent lots, we now use (15) to determine whether there are 
conflicts. If so, we can partition and solve the problem in the following way: 

(i) take the machine which is most conflicted (that is, has conflicts for the largest percentage of time), 

and add it to the partition set. 

(ii) allocate the budget to the parts exactly as was the allocation to the machines belonging to each part 

before. Then, solve the problem for each sub-problem separately. 

Alternately, when we check a new partition we may simply scrap it if there are any conflicts. In 
particular, if conflicts occur when we try to solve without partition, this would imply that we must 
partition the problem, one way or another. To do that, we can use some of the following ways: 

• Partition at machines which conflict while applying the heuristic. (These will probably be slow 
machines.) 

• Try MINPARTIT or medpartit 

• If there are few machines, try all possible partitions. 

To allocate the budget to the parts, we have two alternatives: (i) if we choose to determine the 
partition by the first method, we can allocate the budget accordingly, as discussed above; (ii) in all cases, 
including the former, we can also use dynamic programming similarly to the procedure above. Of course, 
the latter will be more time-consuming. 

Note that if we partition by dynamic programming, allocating the budget to the parts can be done 
together with the partitioning. Otherwise, we start with a given partition, and optimize the allocation of 
the budget to the parts. To do that observe that the savings in the various parts are additive, so a 
knapsack DP algorithm can be used. For this purpose, we have to build table functions for each part 
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showing how much can be saved within it for any given budget allocation. We do this by allocating the 
budget as evenly as possible within each part. If we get conflicts within a part under a certain number 
of parent lots, we look for a different number of parent lots. 

The complexity of the heuristic depends on the choices we make. For instance, if we choose to 
check out all the possible partitions, we have exponential complexity. All the other choices yield 
polynomial complexities. 



9. Dealing with Setups 

(16), below, is a set of conditions which ensure that the setups cannot become critical. If (16) is violated 
for some k, then SUj^ is potentially critical. 

SUj, < MaXj,{SUj,.p + 2j.,p_i.,(Tj + TTj)}; all 2 ^ k < n. (16) 

For every potentially critical setup we have to check the solution for conflicts. If there are 
conflicts, the setups involved are critical indeed. We now list several treatments we can specify if some 
setups are critical. 

The simplest solution is to delay the whole schedule by maXj{D.}, where D- is the amount by 
which SUj conflicts with the original schedule. That is, the most critical setup determines the delay. If 
we do this, however, it may become possible to specify less transfers upstream of the most critical setup, 
and more setups downstream. This may recapture some of the delay within the same budget. 

To minimize the makespan without a budget constraint, we can always partition the problem at 
the most critical machine, and look for a solution which starts feeding it as soon as the setup is 
completed. This may require further partition due to other setups. 

Another approach is to increase Lj (or Yj j or both) enough to ensure that Lj will not reach 
the most critical machine before the setup is finished. 

All the methods above are not guaranteed to produce the minimal makespan. To do that we have 
to resort to the ILP model (see Appendbc). The ILP model itself, however, can be approximated by an 
LP model, thus yielding a solution with polynomial complexity. If we use an LP model, we can obtain 
good integral solutions by the methods listed in Section 10. 



22 



10. Obtaining Locally Optimal Integral Solutions 



To complete the model we still need to adapt the solution of the relaxed problem to the integrality 
constraints. We utilize two procedures: (i) The Feasibility Procedure, which starts with a given make- 
span and produces a feasible integral solution with that makespan--if one exists; and (ii) The Optimizing 
Procedure, which starts with the makespan of the relaxed problem and increases it sequentially until The 
Feasibility Procedure indicates success. The increments by which The Optimizing Procedure incretises 
the makespan are designed to make sure that the minimal feasible makespan will never be exceeded. 
The complexity of the combined procedures is polynomial, and they are fast in practice. 

We assume that the partition of the best relaxed solution is enforced. Our purpose in this section 
is to identify the locally optimal integral solution. That is, the optimal integral solution with the same 
operating procedure and number of lots and sub-lots. Since we do not temper with the partition, we can 
apply our procedures to each part separately. For convenience in presentation, we assume that the 
solution is unpartitioned. 

Let MS* denote the locally optimal integral solution (instead of the globally optimal integral 
solution as above). Then MS* - MS*^ = h* ^ 0 is the difference in makespan due to the integrality 
constraints. Our objective is to find h* and a feasible schedule which makes it possible to process the 
batch within MS*^ + h*. We start with the first operating procedure, and later solve for the second one. 



10.1 Solving for the first operating procedure 



To specify the solution we need to introduce The Feasibility and the Optimizing Procedures as applicable 
to the first operating procedure. 



The Feasibility Procedure; 

Input: Any non-negative value, h. 

Output: Upon success, a feasible integral solution, with makespan MS = MS*^ + h. 

Otherwise, indication that MS* > MS*^ + h. 

Iterations: For i = 1 to k, let 



Lj = min{INT[(MS - T„m + Sj.jCT^ - Tj) - m - Sj.j}) (where = 0); 

Si = Sj_j + Lj; 

Fj = FRACnON[(MS - T„m + Si.i(T„ - Tj) - 

Ej = 1 - Fj (Ej is required for The Optimizing Procedure); 



Success indication: S|^ = m. ■ 
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The procedure simply maximizes Lj subject to the constraint that the makespan should be MS. 
Theorem 1, i.e., specifying that Mj emd should operate continuously, serves to maximize Lj here. 
Note that Fj and E. are similar to fj and bj in Theorem 3, but they are not identical. 

The Optimizing Procedure; 

1. Let hg = 0, j = 1; 

2. iteratively, call The Feasibility Procedure with hj_j; upon success, set 
h* = hj.j and STOP; 

otherwise, set hj = hj ^ + mini{E.}2j,j n.iT.; and start iteration j + 1. ■ 

The basic idea behind The Optimizing Procedure is that if The Feasibility Procedure does not 
produce a feasible solution, there must exist some lot which should include at least one more item. 
Therefore, we look for the smallest addition to h which will cause one lot (at least) to increase by one 
item, and run The Feasibility Procedure again for the new h. We refer to the resulting values of h as 
jump points. A tip to the wise may be in order here: when programming The Optimizing Procedure add 
a small amount, say 10‘^, to h- before calling The Feasibility Procedure. Otherwise, the jump point may 

J 

be missed due to rounding errors. 

The combined procedures’ worst case complexity is polynomial. The proof is a direct extension 
of Theorem 5 in [17]. The procedures were programmed for the two machines case, and the numerical 
experience is that the optimum is usually achieved with less than m/3 iterations. The bound of Theorem 
3 was usually at least twice that of the actual solution. 

10.2 Solving for the second operating procedure 

The complication here is that we have to adjust the sub-lots and the parent lots to be integral. A key 
observation we use to resolve this issue is that processing a parent lot under this scheme is an instance 
of processing a batch under maxpartit. Therefore, if we know how many items are in a lot, we can 
apply the solution of 10.1 to each pair separately. This will provide the locally optimal integral sub-lot 
sizes. 

We proceed to examine how many units can be included in L| for any given makespan. Let /\tj 
measure the time interval between starting Lj on Mj and on then Xj = INT(/\t./2^^jn.jYj j |^T|^) 
is an upper bound on the number of items that Lj can comprise. The Yj j values are obtained as per 
Section 7. When we run The Extended Feasibility Procedure (below), it is possible to compute /\tj for 
each lot. Hence, we can simply try to fit Xj units in Lj, and if necessary decrease it to x* - 1, x. - 2, and 
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so on. The largest feasible value is the one we specify. Finally, let MSj(x) be the minimal makespan for 
a batch of x units under MAXPARTIT with the same number of lots as the number of sub-lots in parent lot 
i between the same machines, and we are ready to state the procedure. 

The Extended Feasibility Procedure; 

Input: Any non-negative value, h. 

Output: Upon success, a feasible integral solution, with makespan MS = MS*^ + h. 

Otherwise, indication that MS* > MS*^ + h. 

Iterations: For i = 1 to k, let 

At; = MS . T„m + S,.,(T„ . Tj) - - SUj; 

L| = max{INT(/\tj/2^,,in.iYj^j^l^T|^), m - 

REPEAT: if MS*(Lj) > /\tj + L-T„ then Lj = Lj - 1; 

UNTIL: Ms!(Lj) < Atj + LjT„; 

Si = Sj_j + Lj (where Sq = 0); 

Ej = MSj(Lj + 1) - A^i ■ (^i required for The Optimizing Procedure); 

Success indication: Sj^ = m. ■ 

The Extended Optimizing Procedure is almost identical to the preceding version. The only 
difference is that when updating h we use hj = hj_j + min{Ej}, i.e., we do not multiply min{Ej} by 

Though slightly more complicated than the algorithm for the first operating procedure, the 
algorithm here is still polynomial and tractable. 

In fact, our method is actually too good in a sense, because in practice we’ll need some time 
buffers (see [14, 17]), to accommodate fluctuations in the processing rate etc. These buffers will probably 
be large relative to the accuracy of the procedure. Therefore, it makes sense to let the increment in h 
be "too large." 



11. Modifications 

So far we assumed that the cost of each transfer is a constant, regardless of the quantity transferred. We 
did not treat the issue of work-in-process explicitely. That is, we tacitly assumed there is enough storage 
space near each machine for any lot size. We did not consider resources that require the same time to 
process any lot size— such as ovens. Finally, we assumed that the product is processed on a single linear 
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sequence of machines. To adapt the model for implementation we may have to deal with some or all 
of these issues. In this section we outline how to do this. 

Suppose the real cost of a transfer from M. is of the form aj + bjL, where a. and bj are positive 
constants and L is the lot size. Than summing for all lots, our cost is a|k + b|m. b«m is fixed for the 
batch, and aj takes the place of our transferring cost, C|. Therefore, all we have to do is reduce the 
budget by m2bj, and our model still applies. Usually we can approximate the real transferring costs by 
such a function, so our model is not restrictive here. 

Next, let us discuss the WIP. There are two issues involved: (i) the money invested in this 
inventory, and (ii) congestion in the plant. The money invested in WIP is only an issue if the raw 
materials and purchased components for some of the items of the batch can be acquired during the 
processing. If this is the case, we should consider using smaller batches. If we insist on large batches, 
however, we can accommodate a restriction on the rate of investment in WIP by specifying a dummy 
machine, Mq, in front of Mj. A transfer from Mq will cost Cq, reflecting the fixed transaction cost of 
ordering/receiving the materials plus the cost of releasing them to production. If the speed assigned to 
Mq is not less than that of the slowest machine, the makespan can still be minimized as before. This may 
require a larger budget, however. 

It is possible to use the model to choose the best speed for Mq so the total cost of WIP and 
transferring will be minimized for any feasible makespan. Next, if we have the value of a savings of a 
time unit in the makespan, we can minimize the total WIP, transferring, and makespan cost. This will 
require using the model as input for a search procedure that will search for the optimal Tq. Note that 
if Tq is set equal to Tp = Max{T.}, then at least between Mq and Mp the model will call for equal 
parent lots. The sub-lots are still likely to vary, however. 

If we can sell the first items of the batch before finishing the processing, than it makes sense to 
allow Mjj to work intermittently too. We can do this by a symmetric dummy machine, j. Again we 
can optimize Tjj^.j similarly to the case of Tq. 

If the WIP is a problem due to congestion, the formal problem becomes tough mathematically. 
We can handle it in practice by dividing the batch and the budget to (roughly) equal sub-batches, and 
solving for each sub-batch by our model. Since the sub-batches are equal, then for every number of sub- 
batches we have to solve the model once, and check the solution for congestion. As a rule, we should 
divide the batch to the largest possible sub-batches which do not cause congestion. 

Next we discuss special resources such as ovens, which take the same time to process a lot 
regardless of its size. A convenient way to deal with these is to model them as transfers, rather than as 
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machines. This raises a sub-problem: make sure that the resource will be available for the next lot in 
time (i.e., the lots must not be too small). The sub-problem can be solved by modifying a model 
presented in [17, Section 5], where the number of vehicles is limited. Our special resource acts as such 
a vehicle, connecting the two adjacent machines. [17]’s model can also serve if there are j such resources 
in parallel, or if the resource has a limited capacity. 

When assembling a batch of products, the assembly operation may be fed by more than one line 
of machines. The problem is to coordinate these lines to feed the assembly on time. To solve this 
problem we can model the assembly as the last machine (M^^) for all the lines feeding it. This creates 
an opportunity to optimize the transportation budget allocation to the sequences. If the objective function 
is to minimize the project makespan, this can be done by The Greedy Heuristic. At each stage, only 
transfers which decrease the makespan of the longest sequence are considered. 



12. Conclusion 



We developed fast solution techniques for the single job-several machines transfer lot sizing problem. 
The techniques can be implemented in new or existing MRP packages. They are easy to program, and 
do not require support by additional mathematical programming modules. We showed that by allowing 
the lot sizes to vary, the number of necessary transfers tends to be 0(log m), and that by allowing 
intermediary machines to idle intermittently the makespan can be decreased considerably. We conclude 
the paper with a partial list of open research questions (see [17] for other open questions). 

• Determine the complexity of the problem; i.e., find a solution in P or prove NP-completeness. 

• Develop more heuristics, as long as no efficient polynomial solution is found. For instance, in the 
second operating procedure we may allow intermittent idling within each parent lot, or allow sub- 
parent lots. That is, use the first or the second operating procedure within the second operating 
procedure. 

• Generalize the problem for several jobs in a flow shop environment (see [1] and [17] for 
preliminary results). 

• Investigate the implications for the Job Shop Scheduling Problem (not necessarily a flow shop); 
combined heuristics. 

• Relax the assumption that the production rates are deterministic and known exactly. This issue 
includes the problem of obtaining the best estimators for the true rates of production to minimize 
the expected makespan. (See [17] for some basic sensitivity analysis results which can be extended 
to the present model.) 
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Consider the case where several operations are required on the same machine, calling for 
intermediary setups. 

Introduce multidimensional budget constraints, e.g., manpower and equipment. 
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APPENDIX: 



This appendix lists the formulation of the problem as an ILP model, and supplies proofs and details 



omitted in the main body of the paper, 

ILP Formulation of Problem IP) : Let ty be the time item j (j = 1,2, ...,m) is transferred to 
(i = l,2,...,n-l), and t • is the time item j finishes processing on M . Let yj : = 1 if t; . coincides with 
finishing the processing of item j on M,, and y; ; = 0 if item j is held until at least one additional item 

1 IjJ 

is processed. This leads to the following ILP formulation: 



of the setups. With k = 0, (A2) ensures that all items will be transferred from M*,j prior to being 
processed by Mj; and with k >. 1, it ensures that item j will not be processed before items j-1, j-2, and 
so on. (A3) is lax if yy = 1, i.e., if a transfer follows the processing of item j on Mj immediately; if 
yy = 0, (A3) implies ty ^ hj + i» target function this will be satisfied as an equality. 

Finally, (A4) is our budget constraint. Note that the number of constraints involved is polynomial, but 
rather large, so the ILP approach may not be attractive in practice. 

For a flow shop environment, it is straightforward to generalize this ILP model for several 
consecutive jobs. 

Theorem 2: If the feasibility conditions are satisfied, [(5)] ^ 1 AND [(5)]Q*^'^ ^ 1, then Lj = [(5)]; else, 
if the feasibility conditions are satisfied, then The Adjustment Procedure yields an optimal solution. 

Proof: We first show that (3) and (5) yield the best super-relaxed solution. If Mjj can start upon 

receipt of Lj, the makespan is Lj(Tj + T 2 +...+ T^ j) + mT^^; i.e., the time required for the first lot 
to reach plus the time required by to finish the batch. By symmetry, if is ready in time for 
the last lot, the makespan is Lj.(T 2 + T^ + ... + T^^) + mTj. Clearly (3) and (5) ensure 



s.t. 




(Al) 

(A2) 

(A3) 

(A4) 



'J 

We have about nm /3 constraints (predominantly (A2)’s), and most of them will be lax. (Al) takes care 



Lj(Tj + T 2 + ... + Tjj_j) + mTjj = L^(T 2 + T^ + ... + T^j) + mTj. We proceed to prove by 



Al 



induction that any other lots are not optimal: 

(a) Let k = 2. If Lj > Lj, then the makespan is at least Lj(Tj + T 2 +...+ j) + > 

Lj(Tj + T 2 +...+ Tjj.j) + mTjj. Alternately, if Lj < Lj, then L 2 > L^, and the makespan is at least 
I-2(T2 + T3 + ... + T^) + mTj > I-2(T2 + T3 + ... + T^) + mTj. 

(b) Let k ^ 3, and the induction assumption is that (3) and (5) are optimal for k-1 lots. For k = 3, 
we proved the induction assumption in (a). If Lj^ > Lj, then the proof in (a) holds here too, and the 
makespan will be larger than for L^. If Lj < L^, then by the optimality criterion of Bellman, the best 
we can do downstream is solve the k-1 case for the remaining m-Lj items. This implies Lj = 
Lj(m-Lj)/(m-Lj) > Lj ; for all j = 2, 3,..., k, i.e., Lj^ > L|^, and the proof in (a) holds again. (Note that 
in the case for which we applied the optimality criterion of Bellman, M 2 and the downstream machines 
will have to idle between Lj and L 2 .) 

It remains to show that The Adjustment Procedure, which we use if [(5)] < 1 OR [(5)]Q*^'^ < 1 
(i.e., the super-relaxed solution is not also the regular relaxed solution), preserves the optimality of the 
solution. We concentrate on Part (a) of the procedure. Part (b) being symmetric. cannot start before 
time 2j=i „.jTj, so by feeding it at this time and making sure it can operate continuously until it finishes 
the batch, the makespan will be minimized. Except for Step 4, this is the case here. As for Step 4, it 
may force some of the last items to be transferred from Mj toward one-by-one, as soon as they are 
finished. In addition, if Step 4 reduces L- , to increase L-, then L- < will be released for trandsfer sooner 

J J J 

than scheduled before. All this cannot cause any delay relative to any feasible solution, so it preserves 
the optimality. Finally, the stopping criterion that is incorporated in Step 4 is valid because the original 
lots are non-increasing. ■■■ 

Proofs Regarding Heuristics 1 and 2: 

We now prove that Heuristic 1 yields a feasible solution that increases the relaxed makespan by 
Max{ej}2j^j „.jTj. First note that S|^ = m, so it does not require rounding. Next, if some of the last lots 
contain one unit each, as a result of Step 4 in The Adjustment Procedure, then the corresponding Sj 
values are integers as well, and do not require rounding either. Therefore, rounding any Sj value up 
cannot cause any subsequent lot to be less than 1. Hence, the lot sizes are feasible. Now, if we start 
processing under this new scheme, will have to wait e|2j,j„.jTj for Lj (relative to the relaxed 
solution), or Max{e.}2j,j„.jT. at most. This completes the proof. 

J 
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The proof that Heuristic 2 yields a feasible solution that increases the relaxed makespan by 
Max{fj}Sj^ 2 ,nTj is by symmetry: rounding up for the original problem truncates the symmetric problem; 
^j= 2 ,nTj assumes the role of and f| replaces e*. ■■■ 

Theorem 4 : For any partition {p(o) = l, p(l), p(2), p(r), p(r+l) = n} (feasible or not), 

MMRi „ > 2i=o.rMMRp(i)^p(i + l). 

Proof : Substituting from (1) and (7), we need to show that 
- Max,=,,„{Tj} 

— ^s=0,i^i=p(s),p(s + l)Ti * ^s-0,r^^i=p(s),p(s + l){"ri}' (^) 

There exists an index w (0 ^ w ^ r) such that MaXj,j„{T|} = MaXj^p^^^p^^+j^lTj}}. Therefore the right 
hand side of (A5) can be written as 



^i=l,n"^i ^S = l,r"^p(s) * ^=0,r^^=p(s),p(s + l){T'i} “ 

2i = l,nTi + Sj^i^^Tp^gj - \ =0.w-l^^i=p(s).p(s+l){'^i} 



- jMaxj^ p(s),p(S + l){T|} - MaX;^j„{T|}. 



Subtract 2;,j„Tj - Maxj^jnlTj} from both sides, and it remains to show 



^ = l,w-l'Tp(54. * ^s=0,w-l^^i=p(s),p{s + l){'^i} ^s=w+l.r'^p(s) " ^s=w+l,r^^^=p(s),p(s + l){'^i} — 

But for any 1 ^ s ^ w-1, Tp^^^ - MaXi^p^j^p^j^j^lTj} <_ 0, and similarly for any w+ 1 ^ s ^ 1, 

Tp(g) - MaXj=p(,)p(j^i){Tj} ^ 0, and the theorem follows. ■■■ 

Lemma 1 : Let (i, j) be any pair such that i < p < j, and such that Tp > Max{T., Tj}, then (i, j) is an 

infeasible pair. 

Proof : Define SUMj = 2j^j+,p.jTg and SUM 2 = Sj^pj.jTg, then 

Qj^p = (SUMj + Tp)/(SUMj + Tj) and Q.j = (SUMj + SUM2 + Tj)/(SUMj + SUM2 + T-). We 
proceed to check if (i, j) is a feasible pair, and clearly if so then we must have Qjj ^ Qj p. After some 
algebra we get Qjj ^ Qj p < = = = > Tj(SUMj + SUM2 + Tj) ^ Tp(SUMj + SUM2 + Tj), but 
Tj < Tp = = = > Tj(SUMj + SUM2) < Tp(SUMj + SUM2), and Tj < Tp = = = > TjTj < TpTj. 
Hence Qj „ > Q; i> ^^d (i, j) is infeasible. ■■■ 

Theorem 5 : Let MINPARTIT= {p(o) = l, p(l), p(2), ..., p(r), p(r+l)=n}, then 

MMRi n - ^i=o,r^^^p(i)^p(i + l)- 

Proof : By observing the proof of Theorem 4, MMRj = 2s=o,ri^i^^^p(s)p(s + l) 

'^pfi+l) ~ ^^=p(i).p(i+i)^^s^ ’ i “ 0> 1> ••••» w-1, AND Tp^j^ — MaXj.p^-j^P^j+j^lTg}, 

for all i = w+ 1, w+2, ...., r, where p(w) is the index of the slowest machine. Now, if the only partition 
is at the slowest machine, the theorem is satisfied trivially. By Lemma 1 Mp(w) is part of the partition, 
so any other machine in the partition must either precede it or follow it. We concentrate on the former. 
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and look at Mp(j) for some 0 ^ i < w. We have to show >. Tp^jj, and then by Lemma 1 we’II 

have shown that Tp^.^.j^ = Max^^p^j^pp+ijITg} as required. If i+1 = w, then this is clearly true, so we 
assume i < w-1. But by construction of the minimal partition we have Qp(j)^p(i + i) > ^p(i) p(w) — ^ 
(since Tp(^) ^Tp(.p = = = > > Tp^p. 

This completes the proof for 1 ^ i < w. As for w < i ^ r + 1, this side follows by symmetry. ■■■ 

In order to prove Theorem 6, we need two additional lemmas, not listed in the main body of the 

paper. 

Lemma Al : Let 1 ^ p(s-l) < p(s) < j ^ n, where (p(s-l), p(s)) is a pair in MINPARTIT, then 
^P(s)J ^p(s-l),p(s)- 

Proof : Let SUMj = 2i=p(j.j)+jp(5).iTj, and SUM2 = 2j=p(5j+,j.jTj, then the lemma states: 

(SUMj + Tj)/(SUM 2 + < (SUM, + Tp,^p/(SUM, + 

We proceed by negation, i.e., assume 

(SUMj + Tj)/(SUM 2 + Tp(j)) > (SUM, + Tp,^))/(SUM, + Tp,^.,,) = = = > 

SUM,Tj + SUMjTp,^,, + TjTp,,.„ > SUM,Tp,^, + SUM^Tp,,, + (Tp,„)2. (A6) 

Now look at Qp(s.i)j - Qp(s-i) p(s)' definition of MINPARTIT, it must be strictly negative. That 

is 

(SUM, + Tp,^) + SUM^ + Tj)/(SUM, + Tp,^, + SUM, + Tp,^.,,) - 

(SUM, + Tp(^p/(SUM, -f Tp(j.,)) < 0 = = = > 

(SUM, + Tp(j, t SUM, + TjKSUM, + Tp,^.,,) ■ 

(SUM, + Tp,j, + SUM, Tp,,.„)(SUM, Tp,^,) < 0. 

By simple algebra, this reduces to 

SUM,Tj + SUM,Tp(^.„ + TjTp,^,) < SUM,Tp,,, + SUM,Tp,,, + (Tp,,,)^ 
thus directly contradicting (A6). ■ ■ ■ 

Corollary Al; Under MINPARTIT, the series {Qp(i) p(i + i))i=o,r monotone decreasing. ■■■ 

Lemma A2 : Under the conditions of Lemma Al, let p(s-l) < i < p(s), then Qj p^^^ > Qp(s)j- 

Proof: By Lemma Al, Qp(g)j < Qp(s-l),p(s)’ suffices to show 

^i,p(s) — ^p(s-l),p(s)- 
By the definition of MINPARTIT, 
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^p(s-l),i — ^p(s-l),p(s)- 



(A8) 



Now look at the symmetric problem, and the symmetric minpartit is the original minpartit, though 
listed in reversed order. This is true because feasibility of pairs is unaltered under symmetry, and 
MINPARTIT is essentially the set of the largest possible feasible pairs. In the symmetric problem, 
^p(s-l),p(s) replaced by VQp(s-i),p(s)’ which applies to the pair (p(s), p(s-l)). Therefore, by 
symmetry to (A8) we obtain 1/Q| <, VQp(s-i) p(s)> which leads directly to (A7). ■■■ 

Theorem 6 : Let PARTIT be any feasible partition which is not identical to minpartit then minpartit must 
be properly contained in PARTIT. 

Proof : Let MINPARTIT = {p(s)}s=o,r+i. and we use simple indices such as i, j for machines in partit 
which are not in minpartit. First, let us show (by negation) that partit cannot be properly contained 
in MINPARTIT. Suppose 1 ^ s ^ r is the smallest index such that p(s) belongs to minpartit but not to 
PARTIT, and let p(t) such that s < t^ r + 1, be the smallest index of a machine which belongs to 
MINPARTIT and to PARTIT (recall p(r+ l)=n, and Mjj is included in all partitions, therefore, such a t must 
exist). By construction of MINPARTIT, Qp^ 5 .j) p^j) < Qp(s-i) p(s)’ hence (p(s-l), p(t)) is not a 
feasible pair. This contradicts the assumption that partit is feasible and contained in minpartit. 

Assume then that PARTIT includes at least one machine which does not belong to MINPARTIT. Pick 
the machine with the smallest index, say i, which belongs to PARTIT but not to MINPARTIT, and clearly i 
> 1 (since Mj = belongs to minpartit); therefore there exists an index s such that 1 < s ^ r+ 1, 

p(s-l) < i, and p(s) > i. Let p(k) be the index of the machine paired with i from below, i.e., the pair 
(p(k), i) is in partit. 

We now show that p(s-l) is p(k). By construction of MINPARTIT, if p(k+ 1) ^ p(s-l) then 
^p(k) i ^p(k) p(k+l) (p(J^)> i) is an infeasible pair, contradicting the feasibility of PARTIT. 

Hence p(s-l) must be in partit. By symmetry, it is clear that if j is the last machine in 
PARTIT - MINPARTIT, then all the machines with a larger index in minpartit must also be in partit. It 
remains to show that no intermediate machines in MINPARTIT are strictly within feasible pairs of machines 
in PARTIT (where each pair may include up to one machine which is also in minpartit; we dealt with the 
case where both of them are in minpartit above by showing that PARTIT cannot be properly contained 
in minpartit). Now, if any other index, say i’, exists in PARTIT such that i < i’ < p(s), take the largest 
such i’, and rename it as i. Hence, i is now the largest index of a machine in PARTIT - MINPARTIT such 
that p(s-l) < i < p(s). Also, we know that p(s-l) is in PARTIT. We proceed to prove that p(s) must be 
in PARTIT as well. To that end it suffices to show that Qj = Maxj„p(j){Qjj} (i.e., no feasible pair 
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exists with p(s) strictly within it). This result is assured by Lemma A2, Now, look for the next machine 
in PARTIT - MINPARTIT which can agtdn be called j. We just proved that the machine in MINPARTIT 
nearest to i from above must be in partit. By symmetry, the nearest machine in MINPARTIT to j from 
below must be in partit also, and it follows that any machines in MINPARTIT between those two are also 
in. Now, rename j as i and repeat the whole procedure until no machines are found in 
PARTIT - MINPARTIT. ■■■ 

Lemma 2 : Let (i, j) be a pair such that T|^ ^ min{T-, Tj} ; for all i < k < j, then (i, j) is a feasible pair. 

Proof : Let SUM = and assume Qy >. 1, then 

Qj = (SUM + Tj,)/(SUM + Tj) j< 1 (since Tj^ :< T-). = = = > Qj |^ <. Qjj; for all k, and the lemma 
is satisfied. If Q: ; < 1, the proof is by symmetry. ■■■ 

Theorem 7 : The series {MRjj(k)}^, 2,3 m “ monotone decreasing. 

Proof : For convenience, we use Q for Q- ; where there is no risk of confusion. MR. .(k) is simply the 
difference between Lj for k-1 and for k transfers multiplied by 2j^ij.jTg. Assuming K-j > k ^ 2 we 
obtain: 



MRjj(k)/ 25 ^ij.,Tg = Lj|k-1 transfers - Lj|k transfers = 

m/(l + Q + Q2 + ... + 0*^-2) . m/(l + Q + Q2 + ... + Q*^-l) = 

mQ*^- V[(l + Q + + ... + Q*^*2)(1 + Q + + ...+ Q*^'!)]. (A9) 

MR. j(K; j) is bounded from above by the value indicated for it by (A9), so it is enough to show that (A9) 
leads to a monotone decreasing series. Assume now that Q :< 1. Under this assumption the numerator 
in (A9) is monotone non-increasing and the denominator is monotone increasing with k. Hence, 
MR.j(k) < MRjj(k-l) as required. The proof for Q > 1 follows by symmetry, since in the symmetric 
problem Q < 1, but the series {MR. .(k) is identical. ■■■ 

Note that (8) can be developed directly from (A9) by using the geometric sum formula where 
applicable. 



Theorem 8 : MR. .(2) > MMR- J2. 

Proon Assume Tj > Tj, = = > Q = Q.^ = I,..,, < 1 and MMR.^ = (m - 

By (8), MR.^j(2) = m(l - 1/(1 Q ) ) 2^ = ij., T . Since Q ^1, 



2MR.j(2) > > (m - = MMRjj. 

By symmetry, the same result holds if Tj < Tj. ■■■ 
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Partitioning the Problem bv Dynamic Programming 



Let: 

• FB = B - 2Cj, be the total free budget after accounting for the essential transfers (i.e., one transfer 
for each pair (i, i+ 1)). 

• F = the free budget remaining for allocation at any stage; e.g., at the first stage, F = FB. 

• TGj(F, j, k) = the total makespan reduction from Mj and downstream if we indicate k transfer 
lots from Mj to Mj, and use the rest of the free budget from Mj and downstream optimally. (In 
this definition and the following ones we assume (i, j) is a feasible pair.) 

• TGj(F) = the total makespan reduction possible from Mj and downstream if we have a free 
budget of F at Mj. 

• "TRj i(k) = the total makespan reduction accumulated for pair (i, j) using k transfers between i and 
j, then, 

TRi .*(1) = 0 (the first transfer is essential), 

TRjj(k) = TRjj(k-l) + MRjj(k) ;2<k<K.., 

TRy(Kij) = MMRy 

• LABELlj(F) = the number of transfers required from Mj when the free budget remaining there 

is F, to achieve the optimal makespan reduction indicated by TG j(F). 

• LABEL2j(F) = the index of the machine paired directly to Mj when we have a free budget of F 

there, to realize the optimal makespan reduction indicated by TG *(F) from Mj and downstream. 

We are now ready to state our recursion formulae. First we have 

TG;.,(F) = TR „,1 „{Min{INT(F/C„.,), K„., „) ). (AlO) 

Next, for stage i, assume we have TGj(F) from the former recursions for all j > i (for i = n-2 we 
have (AlO)), then for any j such that (i, j) is a feasible pair we have 

TGj(F, j, 1) = Tg](F) ; for all 0 ^ F ^ FB, 

TGj(F, j, k) = TRj^k) + Tg!(F - (k - ; any F, 2<k<Kj^, 

Tg!(F) = MaXj,JTGj(F,j, k)}. 
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And thus we can compute all the TGj(F) values, including TGj(FB). Finally, to trace the optimal 
allocation from Mj and downstream, we compute the LABELlj(F) and LABEL2-(F) values during the 
recursion by 

(LABELlj(F), LABEL2j(F)) = args{TGj(F, LABELl.(F), LABEL2j(F)) = TgJ(F)}. 



The Complexity of the DP Partitioning Al gorithm 



Assume we use tables where the row k corresponds to utilizing k transfers in the next stage, and each 

column corresponds to a free budget value. Therefore each table has 0(F)0(max # of transfers) entries. 

There are at most n(n-l)/2 feasible pairs (at least n- 1 , but this is not important for the worst case 

analysis), which can be identified in 0(n ). Thus we have to build 0 (n ) tables, each based on up to n- 

'1 

1 possible routes to feasible downstream pair-mate machines. This leads us to 0(n )0(F)0(max # of 
transfers). Generally 0(F) depends on the budget, but cannot exceed 0(m). As for the max # of 



'I 2 

transfers, if Qj • is likely to be 1, we have to consider up to m transfers, leading to 0(n m ). If Q. . is 
not likely to be 1 more than a bounded number of times which is not dependent on m, than by ( 11 ) and 
(12) we know that K.j is O(log m), leading to O(n^mlog m). 

Note that by using Theorem 6 , we can save a lot of effort when looking for all possible feasible 
pairs that include Mj. The theorem allows us to confine such searches to the subset to which i belongs 
in MINPARTIT. Furthermore, we can partition the problem to the parts implied by MINPARTIT, and use 
a similar master program to assign the budget to the parts. The advantage is realized if FB is larger than 
the budget which can be utilized in some single parts. For instance, if Mp(s) ^p(s+l) 
MINPARTIT, the optimal solution cannot specify spending more in this part than the cost of ^^( 5 ) p(s + l) 
compound transfers, which may be significantly less than FB. Assigm’ng the budget to the parts can be 
solved as a simple instance of the dynamic programming knapsack model (e.g., see [4]). In fact our 
algorithm above is a direct extension of this classic model. 
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Calculating L» and Y - for the second operating procedure 

We now discuss in more detail how to calculate L- and Y- • i^. We follow the schematic outline of Section 
7 step by step. We repeat the outline below, for convenience. 

(1) At stage i, given L| and Y . for p < i and all j, k, calculate Y. . ^ for all j, k; if i = K, go to (3); 
else, set i = i + 1 and go to (2); 

(2) Given L* and Y_ . ^ for p < i and all j, k, calculate l| and return to (1) 

(3) For i = 1 to K let Lj = mL|/2L* 

Step 1: By (13) and (14) we know exactly when L| starts on Mj (STj j), and when its first sub-lot is due 
at Mjj (STj jj). Then the problem of finding the values for Yjj can be solved by applying the super- 
relaxed solution of the two machines model recursively. That is, if parent lot i is sub-divided to k sub- 
lots when processed by Mj, we use (5) and (3) with Q = Tj^.j/Tj and m = 1 (since the Y values sum 
to 1, rather than to m). This policy, if feasible, is optimal when the number of sub-lots for each machine 
is given. The optimality follows directly from the optimality of the two machines model solution. Any 
other choice of the Y- . i, values will lead to unnecessary delays. 

Step 2; Observing the solution for the Y values as discussed in the preceding paragraph, we note that 
they are invariant with Lj. In contrast, the time elapsed between STj j and the instant the first sub-lot 
of Lj can reach is a function of Lj, namely Lj\^, „,j(Yj j j^Tj^). Since we have STj we know the 
time alloted for this purpose, which is [(14)] - [(13)]. On the one hand, if the value we get by the 
tentative Lj is less than the alloted time, then we could have processed more items for the same transfer 
costs. Hence, Lj should be larger. On the other hand, by Theorem 1 we know that M^j does not have 
to wait. Hence Lj must be exactly large enough to use up all the time alloted to it by [(14)] - [(13)]. 
That is 



L, = (STi_„ - STi_i)/J,.,,,(Yy 

Step 3: At the end of Step 2, we have the optimal solution for a batch of 2Lp items. Step 3 simply 
adjusts all the l| values so that their sum will be m. 
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Figure 1. m 250; Tj - 1; T 2 - 2; = 1; = 3. Two equal lots. No idling. 
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Figure 2: m = 250; T. = 1; T 2 = 2; T-, = 1; = 3. Two equal lots. Idling allowed. 





1; = 3. Two varying lots. Idling allowed. 





Figure 5: m = 5; T, = 1; T, = 1; T, = 5; C, = 1; C, = 10; B 






Figure 6: m = 250; T, = 1; T, = 1; T. = 2; T. = 3. Two lots. 




Figure 8: m = 250; T, = l; = 2; = l; = 3. L. = 100 = = > conflicts. 
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